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We study the conductance of normal-superconducting quantum dots with strong spin-orbit scat- 
tering, coupled to a source reservoir using a single-mode spin-filtering quantum point contact. The 
choice of the system is guided by the aim to study triplet Andreev refiection without relying on half 
metallic materials with specific interface properties. Focusing on the zero temperature, zero-bias 
regime, we show how dephasing due to the presence of a voltage probe enables the conductance, 
which vanishes in the quantum limit, to take nonzero values. Concentrating on chaotic quantum 
dots, we obtain the full distribution of the conductance as a function of the dephasing rate. As 
dephasing gradually lifts the conductance from zero, the dependence of the conductance fluctua- 
tions on the dephasing rate is nonmonotonic. This is in contrast to chaotic quantum dots in usual 
transport situations, where dephasing monotonically suppresses the conductance fluctuations. 

PACS numbers: 74.45.-l-c, 85.75.-d, 73.63.Kv, 03.65.Yz 
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I. INTRODUCTION 

The triplet superconducting proximity efFect^'^'^ in 
half-metals (fully polarized ferromagnets, conducting 
only for one spin direction) has received a considerable 
attention recently, both theoretically'^i^i^ii^i^ii^iii and 
experimentallyJ^ii^iii^ii^ The mechanism behind the ef- 
fect is the process of triplet Andreev reflection at the 
half-metal-superconductor interface.^'"* The key ingredi- 
ent that allows and influences this reflection process is 
provided by the magnetic properties of the interface be- 
tween the half-metal and the superconductor: it should 
have a magnetization that is misaligned from that of the 
half-metali^i^ The role of such an interface is to break all 
the symmetries in spin-space, thereby allowing for the 
spin rotations necessary for the triplet Andreev reflec- 
tion. The properties of the interface, however, are not 
easy to manipulate in experiments, which is the reason 
why only a low proportion of half-metal-superconductor 
samples show behavior consistent with triplet Andreev 
reflectionJ^iiiii^ Here we study triplet Andreev reflec- 
tion in a setup that is free of this difficulty. The setup 
consists of an Andreev quantum dot^^ (i-e., a quantum 
dot in contact with a superconductor), coupled to a nor- 
mal, source reservoir via a single-mode quantum point 
contact (QPC), see Fig. [TJ Spin-orbit scattering in the 
quantum dot is assumed to be strong enough that the 
direction of the spin is randomized in much shorter time 
than the typical time i^w of the escape from the dot. 
This allows the dot to effectively play the role of the 
interface. The role of the half metal is played by the 
QPC, which is set to the spin-selective /h conductance 
plateau using a parallel magnetic fieldj^i (For simplicity, 
we assume that the Andreev conductance of the contact 
to the superconductor is much larger than e^//i, which 
makes the transport properties insensitive to the details 
of this contact.) 

A surprising feature of triplet Andreev reflection is that 
despite the randomized spin in the quantum dot, the con- 
ductance of such a fully phase coherent, single-channel 




Figure 1: Sketch of the setup studied in this paper. A nor- 
mal conducting chaotic quantum dot (d) coupled to a su- 
perconductor (S) via a many-mode contact, and to a normal 
source reservoir (N), held at an inflnitesimal voltage V , via a 
single-mode, spin-flltering QPC. Dephasing is introduced by 
coupling the dot to a voltage probe (Vp) via a contact sup- 
porting modes with a tunnel barrier (black rectangle) of 
transparency F,^ per mode. 



system vanishes in the zero temperature, linear response 
regime.^- While current can be passed through the sys- 
tem using finite voltages or temperature, it is natural to 
ask, whether there is still a possibility for transport in the 
zero temperature, linear response limit. In this paper we 
show that there is: relaxing the condition of full phase 
coherence enables the zero-bias triplet Andreev conduc- 
tance to take nonzero values. In the remaining sections, 
our goal is to demonstrate this statement by studying 
the behavior of the triplet Andreev conductance in the 
presence of dephasing in detail. 



II. VOLTAGE PROBE AS A SOURCE OF 
DEPHASING 

We introduce dephasing by coupling the quantum 
dot to an additional normal reservoir, which acts as 
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a voltage probej -'^^'^^i^^'^^i^^ A voltage probe draws no 
net current, but it absorbs and reinjects quasiparti- 
cles without a phase relationship, thereby destroying 
phase coherence. Dephasing due to a voltage probe in 
normal-superconductiiiK structures at zero temperature 
was studied in Refs. [23l l 24f25l . for systems where no triplet 
Andreev reflection could occur. The contact to the volt- 
age probe is characterized by the number of modes (in- 
cluding the spin degrees of freedom) A^^ and the tunnel 
probability per mode F^, which together determine the 
dephasing rate as — N^f^T^S/h, where S denotes the 
mean level spacing of the quantum dot. We consider two 
limits, a voltage probe with a single mode, spin-filtering 
contact, A"^ = 1, and a voltage probe with macroscopic 
number of modes, ^ 1. In the first case the dephas- 
ing rate is controlled by the tunnel probability. In the 
second case, for ^ 1, which will turn out to be the 
regime where the conductance behaves nontrivially, the 
transport properties depend on A'^ and only through 
their product, i.e., the dephasing ratci^ The two lim- 
its considered here represent two types of voltage probes 
that appear in the context of spin dependent quantum 
transport. The probe A^ = 1 is a spin-conserving 
probe, while without further specification, the A'^ 3> 1 
can be either a spin-conserving or a spin-nonconserving 
voltage probe. For the systems studied in this paper, 
there is no need for further specification, as the type of 
the voltage probe is unimportant due to the strong spin- 
orbit scattering in the quantum dot. 

We formulate our problem within the framework of 
the scattering matrix approach. The transport quanti- 
ties of interest are expressed in terms of the scattering 
matrix S at the Fermi energy (the chemical potential 
of the superconductor), relating incoming and outgoing 
modes (including the electron-hole degrees of freedom) 
in the contacts to the normal reservoirs. The currents at 
the contact to the source (s) and the voltage probe {(j>) 
are given by^li^ 



Q/3 



T^y - Tr 



(lb) 



where a, (3 — s,(f), and the index e, h refers to electron and 
hole modes, respectively. The voltages Vp are measured 
from the chemical potential of the superconductor, which 
is assumed to be grounded. The voltage is determined 
by demanding that no current is drawn to the voltage 
probe, Icf, = 0. The conductance, defined by G = Vs/Is 
is given by 



4g = i 



he 



7^5 



, (2) 



where we substituted A^^ = 1- Equation ([2]) is the start- 
ing point for our calculations. In what follows, we con- 
centrate on systems where the motion inside the quan- 
tum dot is chaotic. We are interested in the statistics of 



the conductance, which we obtain using Random Matrix 
Theory^^ for the scattering matrix S. 



III. DEPHASING DUE TO A SINGLE MODE 
VOLTAGE PROBE 

We first discuss the case voltage probe with A^^ = 1. 
The calculational advantage of this case is that it allows 
for a problem with minimal dimension, with the single 
mode source contact and a single mode voltage probe 
contact resulting in a 4 x 4 scattering matrix. The par- 
allel magnetic field together with the strong spin-orbit 
scattering places the quantum dot in the unitary symme- 
try class.'^'?^ Consequently, the dot-superconductor system 
belongs to class D in the symmetry classification of Alt- 
land and Zirnbauerf^ which translates to S* = SiS**!]! as 
the only constraint for the scattering matrix, besides uni- 
tarity. (Ej denotes the j-th Pauli matrix in electron-hole 
space.) Assuming that the contact to the source reservoir 
is ideal, the two single-mode QPC-s can be characterized 
by the reflection matrix 















(3) 



where the block structure reflects a grading according to 
the normal contacts and ^ is the reflection phase shift 
for electrons at the voltage probe contact. The statis- 
tical properties of the conductance follow from the dis- 
tribution of the scattering matrix, which is given by the 
generalization of the Poisson kernel^^' 



P{S) oc |det(l-rtS')|^ 



(4) 



The probability distribution is understood with respect 
to the invariant measure dfi{S) on the manifold A4d de- 
flned by = SiS"*!)! in the space of 4 x 4 matrices. We 
parametrize S as 



T12 

Tt e' 



W 
W 



(5) 



where T e (0,1), ^1,^2 S (0,27r), W € SU(2) and 
T — 1(72 • (fj denotes the j-the Pauli matrix in spin- 
space.) The matrix structure in ^ corresponds to 
electron-hole grading. The above parametrization can 
be obtained from the polar decomposition introduced in 
the Appendix. Eq. ([5]) implies that det(S') = 1 and that 
the matrix S'^'^ {S'^'^)'^ has a twofold degenerate eigenvalue 
T. This is true for the generic setups with vanishing lin- 
ear conductance in the fully phase coherent limit, i.e., 
if the closed Andreev quantum dot has no energy level 
at the Fermi energy»il (For a detailed discussion of this 
point we refer to the Appendix.) Using the Euler angle 
parametrization for W , 



e-^(v+^}/^ cos(6'/2) -e^('>-'^)/2 sin(6i/2) 
g«(<p-^A)/2 sin(6i/2) e*('^+'^)/2 cos(6'/2) 

iip,,p,9) e [0,27r] X [0,47r] x [0,7r], 



(6) 
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the invariant measure on is given by dfi{S) oc sin(6'), 
and the conductance in units of e^/ft. is 



G(T, 9) =4 



1 



T sin' (0/2) 



(7) 



The distribution of the conductance is given by 
Pr^{G) = J dn{S)P{S)6{G - GiT,e)), which can be re- 
duced to 



2a' 



p3 „4/G-l 

-Pr^(G) = ^y^ '^^^4/G-x)2(a2- 62)5/2' 



i-(i-r^ 
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(8) 
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for < G < 2 and otherwise. A closed form expression 
can be given for = 1, 



For < F^ < 1, we evaluated the integral ([8|) numer- 
ically. The resulting distribution is shown in the top 
panel of Fig. [2] for several values of F^. In the absence 
of dephasing, the conductance vanishes, corresponding to 
Pr^=o{G) = S{G). With the gradual introduction of de- 
phasing, G is enabled to take nonzero values, leading to 
a widening of the conductance distribution with increas- 
ing dephasing rate, eventually reaching the distribution 
([9|) for F0 = 1. In the bottom panel of Fig. [2] we show 
the dependence of the average and the variance of the 
conductance on F^. While the average conductance in- 
creases monotonically with increasing dephasing rate, the 
variance increases to a maximum at F^ « 0.8, which is 
followed by a slight decrease. The finite value of the con- 
ductance fluctuations at F^ = 1 (corresponding to the 
maximal dephasing for iV^ — 1) indicates that a single 
channel voltage probe can not lead to a complete loss of 
phase coherence - as we will see below, without phase 
coherence, the conductance fluctuations are suppressed 
back to zero. For weak dephasing, F^ <C 1, the conduc- 
tance distribution is rapidly decaying away from G — 
and it has the scaling form Pr^(G) = /(G/F0)/G. This 
results in the dependence 



(G")cxF^, F^«l 
for the n-the moment of the conductance. 



(10) 




GK/h] 




Figure 2; Upper panel; Probability distribution [Eq. ([8}] of 
the conductance for A^^ = 1, for various values of the de- 
phasing rate g^, = T^S/h. The curves, in order of decreasing 
maximum, correspond to F,^ = 0.2, 0.4, 0.6, 1, respectively. 
The empty squares represent smoothed histograms obtained 
from 3000 numerically generated scattering matrices for each 
value of r^, for a system where the superconducting contact 
supports 50 propagating modes. Lower panel: the average 
(solid line) and the standard deviation (dashed line) of the 
conductance as a function of g^. The crosses are results of 
the numerical simulation. 



be possible to make some analytical progress using the 
Poisson kernel distribution of Ref. [32 an d following sim- 
ilar steps to the calculation of Ref. [22, we resort to a 
simpler approach and obtain the statistics of the con- 
ductance by generating an ensemble of scattering matri- 
ces numerically. The scattering matrix S is expressed in 
terms of the electron scattering matrix 



S 



N 



r t' 
t r' 



(11) 



IV. DEPHASING DUE TO A MULTIMODE 
VOLTAGE PROBE 

Now we turn to the case of the voltage probe with 
macroscopic number of channels, TV^ ^ 1. While it might 



of the normal region. Here r is describes reflection from 
the dot through the normal contacts, r' describes reflec- 
tion through the superconducting contacts, t corresponds 
to transmission to the superconducting, and t' to the nor- 
mal contacts. The necessary blocks of S in electron-hole 
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grading are given b y^^i'^^ 

S^^ = r - tV2r'V2(l + rVar'Va)-^ , (12a) 
= -t'V2(l + r'a2r'*a2)-H , (12b) 

The scattering matrix Sn can be expressed using the 
statistical mapping^i^^ 



Sn = VT~^-Vr 



1 



SoVf, 



(13) 



where So is unitary and F is a diagonal matrix contain- 
ing the transmission probabilities of the contacts with 
Fii = 1 corresponding to perfect transmission through 
the single mode QPC and Tjj = for 1 < j < A^^+l de- 
scribing tunneling at the voltage probe. We took Tjj = 1 
for j > + 1, corresponding to the contact to the super- 
conductor. The results do not depend on this choice, as 
long as the Andreev conductance of the contact is much 
greater than e^/h. Using the mapping the distri- 

bution of Sn is obtained by taking 5*0 from the circular 
unitary ensemble^^i^ which we generated numericallyi^ 
For a mutual test of the program and the calculations, 
we first show results for = 1 in Fig. [51 As it is seen, 
the agreement between the calculation and the numer- 
ics is perfect. The conductance distribution in the limit 
^ 1 is shown in the top panel of Fig.[3]for several val- 
ues of the dephasing rate g^. The distribution Pg^ (G) ini- 
tially widens from Pg^=o = <5(G) with increasing and 
then it gradually narrows again to Pg^=oo — ^{G~G class), 
where 



1 



1 



G, 



QPC 



G 



NS 



G, 



QPC 



= 1 



(14) 



is conductance of the single mode QPC and the Andreev 
conductance of the superconducting contact in series, in 
units of e^/h. The dependence of the average and the 
variance of the conductance on g^ is shown in bottom 
panel of Fig. [31 While the average conductance increases 
monotonically to its classical value, the conductance fluc- 
tuations display nonmonotonic behavior, corresponding 
to the initial widening and the final re-narrowing of the 
conductance distribution. Fig. [31 also shows a compar- 
ison between the limits = 1 and ^ 1. For a 
given value of g^, the conductance distribution close to 
G = is suppressed for ^ 1, in contrast to the single 
channel case, where P(G — 0) is finite. The average con- 
ductance increases faster for A^^ ^ 1 towards its classical 
value, while the conductance fluctuations are suppressed 
compared to Na, = 1. 



V. CONCLUSIONS 

In summary, we have studied in detail how dephas- 
ing due to a voltage probe enables a nonzero value for 
the zero temperature, zero-bias triplet Andreev conduc- 
tance in Andreev quantum dots with a single-channel 
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GK/h] 




Figure 3: Upper panel: Probability distribution of the con- 
ductance for ^ 1, for different values of the dephasing rate 
g^. The solid curves, in order of increasing position G of the 
maximum, correspond to g^h/5 = 0.05, 0.5, 1.5, 10, respec- 
tively. The curves are obtained by smoothing histograms from 
3000 numerically generated scattering matrices for each value 
of g^, with = 100, for a system where the superconduct- 
ing contact supports 50 propagating modes. For comparison, 
we show the distribution for A'^^ = 1, g4,h/5 = 0.5 (dashed 
line). Lower panel: the average (solid line) and the standard 
deviation (dashed line) of the conductance as a function of 
g^,. For comparison, the dotted lines show the corresponding 
functions for Ns = 1. 



spin-flltering source point contact. We focused on sys- 
tems where the quantum dot is chaotic, and obtained 
the full distribution of the conductance as a function of 
the dephasing rate for two limiting cases for the number 
of modes in the voltage probe contact, = 1, and 
3> 1. Compared to chaotic quantum dots in other 
transport situations, our findings for the conductance are 
quite unusual. Dephasing is known to monotonically sup- 
press the conductance fluctuations, in general i^^i^^'^''i'^^ 
In contrast, as dephasing gradually enables transport, 
the fluctuations of the triplet Andreev conductance are 
initially enhanced, which is followed by a suppression for 
strong dephasing, i.e., the overall dependence on the de- 
phasing rate is nonmonotonic. 

It is worthwhile to point out that in the A^ 3> 1 case. 
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unlike Ref . |37|, we did not intend to use the voltage probe 
to model dephasing processes intrinsic to the quantum 
dot, since accounting for the temperature dependence of 
such processes would necessitate considering the effect of 
thermal smearing.'^* Instead, our results apply to the sit- 
uation where the dephasing rate is controlled by a real, 
physically present voltage probe. Experimental control 
of the dephasing rate using a voltage probe was demon- 
strated very recently in the work of RouUeau et al^ This 
makes us believe that, in principle, it is realistic to test 
our predictions in experiments. 
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Appendix A: ELECTRON-HOLE SYMMETRY, 
POLAR DECOMPOSITION AND ANDREEV 
REFLECTION EIGENVALUES 

The Andreev reflection eigenvalues Tj are the eigen- 
values of the matrix 5''"^ (5''"^)^. We prove here the 
consequences of electron-hole symmetry on these quan- 
tities, and relate them to the condition of the absence of 
energy level of the closed Andreev quantum dot at the 
Fermi energy. 

Theorem: At the Fermi energy, the degeneracy dj of the 
Andreev reflection eigenvalue Tj is even if Tj ( 1 — Tj ) 7^ 0, 
and det(S') = (—1)'*", where du is the degeneracy of the 
unit Andreev reflection eigenvalue, if present, (i„ = oth- 
erwise. Furthermore, the scattering matrix at the Fermi 
energy can be decomposed in electron-hole grading as 



S = 



U 



R pT\ (V 

V"* ' ' 



^0 {77 \pT R 
where U and V are unitary matrices, 

J 
3 

and p — @j Pj, where 



(A2a) 
(A2b) 



Pi = 



ifT,(l-T,) = 
I T otherwise. 



(A3) 



Proof: Following from the electron-hole symmetry 
S = EiS'*Ei, the scattering matrix has the block decom- 
position 



gee (^ghey 
ghe (^geey 

We start with the singular value decomposition 
S^'- = U'RV' 



(A4) 



(A5) 



where U', V are unitary matrices. Using {S^SY'^ — 1 
and {SS^'y = a, one finds that 

ghe ^ u'*zfV'. (A6) 

Here Z is a block diagonal unitary matrix, 

Z = , dimZj = dj . ( A7) 
i 

Substituting (|A5l) and (|A6| into {S'< S)''" = leads to 



T,il-T,)Z,=-JT,il-T,)Z^ 



(A8) 



For Tj(l — Tj) 0, Zj is antisymmetric, and due to 
its unitarity det(Zj) ^ 0, from which it follows that dj 
is even. Being antisymmetric and unitary, Zj can be 
decomposed as^Siil 



(A9) 



where Uj is unitary. For Tj = 0, 1, Eqn. (|A8|) is auto- 
matically satisfied, without further requirements for Zj. 
For the zero Andreev refiection eigenvalue, if present, we 
can set Zj = UjUj with an arbitrary unitary matrix 
Uj, as for Tj = 0, Zj drops out from (|A6p . For the 
unit Andreev reflection eigenvalue, if present, we write 
Zj = UjU'j with J7j, Uj unitary. Taken together, the 
matrix Z can be written as 



Z = ^Ujp,U^, 



(AlO) 



(Al) where C/j = Uj for Tj ^ 1. Writing U' and V as 



U' = U \^U, 

with J7, V unitary, one finds 

S"^ = URV 
ghe ^ u*pfy^ 



(Alia) 



(AUb) 



(Al2a) 
(A12b) 



which gives the decomposition (jAip upon substitution in 
Eq. (|A4p . The decomposition (|Aip satisfies the unitar- 
ity and electron-hole symmetry requirements, therefore 
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there are no further relations between the matrices U 
and V. The result det(5) = (-1)'*" follows straightfor- 
wardly. □ 

Note that in Eq. ([5]), we assumed det(S') = 1, however 
only det(S') = ±1 follows from S = EiS'*Ei. We show 
below that this is a valid assumption in the generic situ- 
ation that there is no energy level of the closed Andreev 
quantum dot at the Fermi ener gy. Using the channel 
coupled model employed in Ref. l3ll . the scattering ma- 
trix can be expressed as^S 

= i±^. (A13) 

Here the Hermitian matrix He = — Siifl^Si is a pro- 
jection of {TC — E)~^, where H models the closed An- 
dreev quantum dot. If Ti. has no zero eigenvalues, i.e., 
there is no level at the Fermi energy, the matrix He 
can be taken at E = without complications. Follow- 



ing from the symmetry of He=0: the eigenvalues of S 
come in complex conjugate pairs, therefore, det(S') — 1. 
(If there is a level at the Fermi energy, an eigenvalue 
of He can diverge as E ^ while an other can tend 
to zero, leading to a (1,-1) eigenvalue pair of S, and 
thereby to det(S') = —1.) This result, together with 
det(S') = (— 1)'^", contains as a special case the find- 
ing of Ref. [ll|, that for a single mode system, 5''"^ = 
at the Fermi level, if the closed Andreev quantum dot 
has no level at the Fermi energy. Indeed, for such a sys- 
tem, S' is a 2 X 2 matrix, i.e., there is a single Andreev 
reflection eigenvalue. As it is singly degenerate, it can 
be only zero or unity, and det(5') = 1 guarantees that 
it is zero. For the 4x4 matrix in Eq. ([5]), the degener- 
acy of the Andreev reflection eigenvalue also follows from 
det(S') = 1. If there was no degeneracy, the eigenvalues 
could be only a zero and a unit eigenvalue. This would 
mean det(S') = —1. 
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